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Abstract 

We explore the character of a novel set of "field-linked" states that were predicted in [A. V. 
Avdeenkov and J. L. Bohn, Phys. Rev. Lett. 90, 043006 (2003)]. These states exist at ultralow 
temperatures in the presence of an electrostatic field, and their properties are strongly dependent 
on the field's strength. We clarify the nature of these quasi-bound states by constructing their wave 
functions and determining their approximate quantum numbers. As the properties of field-linked 
states are strongly defined by anisotropic dipolar and Stark interactions, we construct adiabatic 
surfaces as functions of both the intermolecular distance and the angle that the intermolecular axis 
makes with the electric field. Within an adiabatic approximation we solve the 2-D Schrodinger 
equation to find bound states, whose energies correlate well with resonance features found in fully- 
converged multichannel scattering calculations. 
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I. INTRODUCTION 



In the modern world of physics, manipulation of quantum phenomena in atoms and 
molecules forms the basis for future applications. With the development of new techniques 



for cooling and trapping polar molecules, new opportunities to harness them appeared 

. In particular, the interactions between pairs of molecules are likely to be 
susceptible to manipulation in an electric field. This in turn may imply an ability to direct 
the course of chemical reactions |ioj, to influence the many-body physics of degenerate Bose 
or Fermi gases composed of polar molecules l], 12, [l^ 14,|l5|, or to manipulate quantum 

bits [lel. 

A particularly attractive opportunity for controlling intermolecular interactions emerges 
in a set of novel long-range bound states of molecular pairs In the presence of an 

external electric field, the counterplay between Stark and dipole-dipole interactions generates 
shallow potentials that are predicted to support bound states of two polar molecules. For 
OH molecules we have estimated that the bound states do not exist at all for fields below 
about 1000 V/cm 17[. Thus the field plays an essential role in binding the molecules into an 
[OH] 2 dimer; we have accordingly dubbed this new kind of molecular state a "field- linked" 
state. The purpose of this communication is to further clarify the structure of field-linked 
(FL) states. Interestingly, quadrupolar interactions between metastable alkaline-earth atoms 
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exhibit similar states in the presence of magnetic fields 

Schematically, the FL states originate in avoided crossings between a pair of potential 
energy curves: one that represents an attractive dipolar interaction converging to a high- 
energy Stark threshold; and one that represents a repulsive dipolar interaction converging 
to a lower-energy threshold. The characteristic size of the FL states is therefore roughly 
determined by equating the dipolar energy fi'^/R^ to the Stark energy fiS. Here R is the 
distance between the molecules, fi is their dipole moment, and S is the field strength. The 
length scale of the avoided crossing is then Rscaie = {ft/^Y^^^^ ~ 1250^^^/3 a "typical" 
dipole moment of 1 Debye, where Rscaie is measured in units of oq (the Bohr radius) and S 
is measured in V/cm. Thus for a reasonable-sized laboratory field of 10'^ V/cm, the size of 
the FL state is ~ 60ao, although extremely weakly bound states can be far larger than this. 

Ref. jl?! described the FL states in this simple curve-crossing picture. Adiabatic potential 
curves for the OH-OH interaction were constructed by expanding the relevant potential 
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into partial waves in the intermolecular coordinate. For clarity, only the lowest partial 
waves, L — 0,2 were included. While intuitively appealing, this picture is inadequate, and 
indeed a partial wave expansion is inappropriate, for the following reason. The dipole-dipole 
interaction can strongly couple different values of L, with a strength on the order of ~ /i^/ R^. 
At the typical scale distance Rgcaie, the dipole coupling exceeds the centrifugal interaction by 
a ratio 2inijL^ /h^Rgcaie, where m is the reduced mass of the molecular pair. For our example 
case oi fjL — 1 D, £ — 10^ V/cm, and for a hght molecule (hke OH) with a reduced mass 
m — 10, this ratio is already 100. The ratio becomes even larger in a stronger field, or 
for a heavier molecule. Therefore L is no longer a good quantum number for the FL states, 
but rather the relative orientation of the molecules is of more significance. 

Accordingly, in this paper wc present a formulation of FL states in terms of potential 
energy surfaces in {R,9), where 9 is the angle that the intermolecular axis makes with re- 
spect to the electric field. Within an adiabatic representation, we compute FL states as 
bound states of a single surface. Quahtatively, these identify the FL states as confined to a 
narrow range about ^ = 0, so that their motion consists primarily of vibration along the field 
axis. Additionally, we show that the binding energies predicted by this adiabatic approx- 
imation agree remarkably well with resonance positions determined from fully-converged 
multichannel scattering calculations. 

II. MODEL 

Because the FL states are generated primarily by the competition between Stark and 
dipolar interactions, our model will focus almost exclusively on these two terms in the 
Hamiltonian. In particular, our simplifying assumptions here are: 

1) The individual molecules are assumed to be in their electronic ground states, to be 
rigid rotors, and to lie in their rotational ground states. It is assumed that none of these 
degrees of freedom can be excited at the large intermolecular separations and low relative 
energies that we consider. 

2) Each molecule is assumed to have total spin j and to have a non-E electronic ground 
state that can support a lamb da- doublet. Again, at the intermolecular separations, energies, 
and fields of interest, it is assumed that j is approximately conserved. We ignore hyperfine 
structure in the model, so that j is an integer for bosonic molecules, and a half-integer for 
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fermionic molecules. While hyperfine structure is well-known to be important in ultracold 
collisions, it is not germane to the main discussion of dipolar interactions, and can in any 
event be included in a straightforward way later. 

3) The projection of each molecule's angular momentum onto its own interatomic axis, 
denoted uo, takes only the two values ±j. As a point of comparison, the energy difference 
between the j = 3/2, Ic^l = 3/2 ground state and the j = 3/2, = 1/2 excited state of OH 
is 270 K I22], so this restriction is not such a bad one. 

4) We work in the limit of large electric field, i.e., in the linear Stark regime where the 
electric field interaction dominates the lambda-doublet splitting. Thus the molecular states 
are characterized by the signed quantities ±uj, rather than linear combinations of u and —uj 
characteristic of the zero-field limit. We will describe some effects of lambda-doubling in the 
following, but they will be perturbative in this limit. A readable account of molecular wave 
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functions in this approximation is given in Ref. 

5) Finally, we assume that the molecules never get close enough together for short-range 
interactions such as hydrogen bonding, exchange, or chemical reactions, to contribute. In 
addition, we neglect long range interactions such as dispersion and quadrupole-quadrupole 
interactions, as being negligible compared to dipole-dipole interactions. 

Although this model does not describe any particular molecule, it lays the groundwork 
for constructing FL states for any desired molecule. To keep the magnitudes of observable 
quantities realistic in the following, we use as model parameters the dipole moment (1.68 
D), lambda- doublet splitting (0.055 cm~^), and mass (17 amu) of the OH radical. 



A. Basis set 



Within the simplifications outlined above, the internal state of an individual rigid-rotor 
molecule is specified by three quantum numbers: j, u, and the projection of the molecule's 
angular momentum on an appropriate external axis, to describe the Stark interaction this 
axis is conveniently taken as the electric field axis. However, to describe FL states we choose 
instead to quantize this angular momentum along the intermolecular axis. This emphasizes 
the dimer nature of the FL states and allows a reasonable description of how the dipole-dipole 
forces act ultimately to keep the molecules from crashing into one another. 
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Each molecule {i — 1,2) is thus described by a rigid rotor wave function, 



{iiljkiUJi) = ^^^^^Dl*^.{ai, A, 7i), (1) 

where notation for the electronic wave function is suppressed, under the assumption that 
it plays no role at the temperatures and electric fields of interest. Here ki and coi are 
the projections of total angular momentum j onto the intermolecular axis and onto the 
molecule's own body-frame axis, respectively. The Euler angles = {ai,Pi,ji) are referred 
to the intermolecular axis. We further couple the molecular spins into a total spin J: 

(ei,e2|(l,2)j;^) = 

J2 (ei \jkiuji) (e2 \jk2UJ2) {jikij2k2\JK) . (2) 

kik2 

Here we introduce the shorthand notation (1,2) do denote the internal molecular quantum 
numbers {jiUi, j2UJ2)- 

As for the relative motion of the molecules, we wish to avoid an expansion into partial 
waves, as mentioned in the Introduction. We thus consider a basis set for the complete wave 
function 

*a!2)ji.(^,^,0,ei,e2) = -^exp{tMcl>)F/^2)jKiR:0) 



X (ei,e2|(l,2)JX), (3) 

where the F's are as-yet-unspecified functions of {R, 9). The projection of the total angular 
momentum onto the electric field axis, A4, is the only rigorously conserved quantity in the 
Hamiltonian for FL states; we therefore separate it at the outset. It will affect the functions 
F via centrifugal energies. 

In addition, the wave functions must incorporate the proper symmetry under the exchange 
of identical molecules, denoted by the operator Pi2- The symmetrized states are constructed 
in Appendix A, and define a pair of quantum numbers s and x: 

P TpM,S _ „pX,s 

A2|(l, 2) JiT), = x\{l,2) JK),. (4) 

The quantities s and x are not separately conserved by the Hamiltonian, but must satisfy 
the constraint 

{-1-1 for bosons 
(5) 
—1 for fermions. 
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Finally, it is useful to consider the effect of the parity operator / that inverts all coor- 
dinates through the system's center-of-mass. Eigenvalues e of this operator are obviously 
not conserved by the electric field, yet we can construct basis sets that are eigenfunctions 
of /, as is done in Appendix A. When we consider matrix elements of the electric field and 
dipole-dipole Hamiltonia, we find that the quantity 

q^ts{-l)^ (6) 

is rigorously conserved (see Appendix B). Our completely general basis then takes the form 

^t^,%ii,2)jK = -^eM^M<P)F^^,lj^\{l,2)JK)^,„ (7) 

whose explicit representation in terms of unsymmetrized basis functions is given in Appendix 
A. 



B. Hamiltonian matrix elements 



To uncover the joint motion in {R, 9) that governs the FL states, we will expand the total 
wave function into the basis and integrate over all other degrees of freedom to derive a 
set of coupled-channel differential equations for the functions F. In this section we therefore 
construct the Hamiltonian matrix elements in the "internal" basis |(1, 2) Ji^)^; g. 

Ignoring the exchange, quadrupole-quadrupole and dispersion interactions as we did in 
jl7|. our model Hamiltonian can be written as 

H=Y.{t, + H!) + y,,. (8) 

i=l,2 

where Tj and Hf are the translational kinetic energy Stark energy of each molecule, and V^^ 
is the dipole-dipole interaction. 

In the following subsections we list the matrix elements of the various terms of the 
Hamiltonian in the unsymmetrized basis. Transformation into the symmetrized basis set is 
accomplished in Appendix B. 



1. Stark Interaction 



An electric field with strength £ that points along the positive z-axis in the laboratory 
frame will have spherical components £q in the reference frame that rotates with the inter- 
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molecular axis. The relation between the two is given by a Wigner rotation matrix: 

S, = SD',^i<P,e,0). (9) 

The components of the molecular dipole moment /x can be written in terms of reduced 
spherical harmonics Ciq(a,/3) where, as above, a and (3 are Euler angles relative to the 
intermolecular axis. The Stark Hamiltonian for a single molecule is then 



(10) 



The integration over each molecule's internal coordinates yields, for the unsymmetrized 
basis set (and remembering that ji = j2 = j ) 3, 

{{1,2)JK\H^\{1',2')J'K') = 










J 1 










X 










I 


-uJi 





><^2 + J+J' 




(11) 

(12) 



where [y] = \/2y + 1. 



2. Dipolar Interaction 

The dipole-dipole interaction reduces to a particularly simple form in the rotating frame: 

. . _ Ml • - 3(^- Hi){R- M2) 

= -^[/^i®Ho- (13) 

Here [/ii ® /i2]o is the (2,0) component of the second-rank tensor formed by the product of /ii 
and fi2- The zero refers to the cylindrically symmetric component around the intermolecular 
axis. 
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Following a treatment similar to the Stark effect above, we note that 

[/ii ® Ho = /i^5I'^i'?('^i'/^i)'^i-'?('^2,/32)(igi - g|20). 



(14) 



Now the angular integration over each molecule's internal coordinates is again straightfor- 
ward, yielding 



{{l,2)JK\V,,\{l\2')J'K') 

.2 




K' —UJ1—UJ2 



( 



3 1 3 

\^ —UJ2 UJ2 

J' 2 J ' 
3 1 J ^ • 

^3 1 J , 

This matrix element is independent of the orientation 6', as it must be. 



(15) 



3. Kinetic energy 

The centrifugal Hamiltonian in the rotating frame is no longer diagonal, but rather couples 
states with K,K ±1 projections. Within our basis it is more convenient to present the 



angular momentum operator as 
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P = j^ + j^- 2J, - {J_J+ + J^J_) (16) 



r\j{i,2)K> = j(j+ 1)1 j(i,2)is: >, 



Knowing that 



andJ±|J(l,2)ir > = ^ J{J + 1) - K{K ±1) 

x|J(l,2)ir±l > (17) 

and using the definition of and F± we have 

xe'^'^(^\J{l,2)K > Ao{K) + |J(l,2)i^- 1 > A^i{K) 
+ 1 J(l, 2)K+1> A^.iK)] F/^,.jj,{R, 6), 



where 



92 d 
+J(J+ 1) - 2K'^ 



A^,{K) = -—^J{J+i)-K{KtI) 



x(—^±M-^] (18) 



de sm{e). 

For convenience, we will in the following neglect the Coriolis-type couplings A±i. Like many 
other perturbations, these can be incorporated later, if necessary. 



C. Schrodinger equation 

Within our scheme we have the following Schrodinger equation 

+ + Vi,,,(R,0))h(n.,l>) = (19) 

i' 

where i = {M, (1, 2)JK, x, q} and fi = Fi/R . 

Solutions of the coupled-channel partial differential equations ()19p. subject to scattering 
boundary conditions, yield both the energies and resonance widths of the FL states. To clar- 
ify the nature of these states, however, we first invoke a Born-Oppenheimer approximation. 
Thus we will diagonalize the model Hamiltonian for fixed values of the pair [R, 6) , and seek 
bound state in one of the resulting potentials. In a single adiabatic surface \/<^'^«"^(/^^ 5))^ the 
Schrodinger equation reads 

_^yadiab^-^^ ^) - f {R, 6) = (20) 

III. CHARACTERISTICS OF THE FIELD-LINKED STATES 

For concreteness, we consider here a pair of bosonic molecules with j = 1, and parameters 
corresponding to the OH radical, as discussed above. From the similar model in Ref. jl7|. 
we then expect to see a small number of FL states at modest electric field values. Our aim 
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100 150 200 250 

R(a.u.) 

FIG. 1: A "slice" through the adiabatic potential energy surfaces for an electric field strength 
£ = 10^ V/cm. In these surfaces g = — 1, wi = 0^2, and even (solid line) and odd (dashed line) 
values of (J) are distinguished. 

in this section is to describe these states approximately in terms of the quantum numbers 
in our basis set defined in the previous section. 

A. Adiabatic surfaces 

The number of adiabatic potential surfaces is set by the number of internal states of the 
molecules. (Contrast this to an expansion in partial waves, where the number of channels is, 
in principle, infinite.) For a pair of j = 1 molecules, the present model contains 36 channels, 
hence 36 surfaces. Moreover, conservation of q implies that these 36 surfaces split into two 
sets of 18 channels each. The surfaces for g = 1 and g = — 1 are identical, if only the Stark 
and dipolar interactions are included, as we assume. We find that including the lambda- 
doublet interaction leaves the q = —1 surfaces unchanged, but introduces some weak avoided 
crossings among the g = 1 surfaces. Since lambda-doubling is a perturbation for the fields 
we consider, we will ignore this small effect. Hereafter we report on the g = — 1 surfaces. 
Additionally, the quantities uji are conserved in the absence of lambda-doubling, meaning 
that we can further classify the surfaces according to whether uji = uj2 or Ui = —002- 

Subdividing the surfaces in this way yields nine surfaces with q = —1 and ui = UJ2, 
which are of greatest interest here. Slices through these surfaces at a fixed angle 6 = 5° are 
shown in Figure ^ Here we take the applied electric field strength to be £^ = 10^ V/cm. 
Empirically, we find that surfaces with even values of (J) (solid lines) are only weakly coupled 
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FIG. 2: Adiabatic potential energy surfaces corresponding to the highest-lying slice in Fig. 1. 
Shown are the 7W = (a) and Ai = 1 (b) cases. 

to surfaces with odd values of (J) (dashed lines). This consideration further reduces the 
number of surfaces necessary to describe the FL states. 

The FL states are bound states of the highest-lying surface in Figure ^ which is clearly 
generated by avoided crossings. The complete surface in the {R, 6) plane is shown in Figure 
121 for both the rotationless case M. = Q and a rotating case with M. = 1. Addition of the 
centrifugal energy makes the M. = 1 surface substantially more shallow than the M. = Q 
surface; in fact we find six bound states for = 0, and only two for M. = 1 (see Table I). 

To gain a better understanding of the nature of the FL states, it is useful to evaluate mean 
values of the quantum numbers in our basis set. In general, the symmetry-type quantum 
numbers x, s, and e are badly nonconserved, and average to zero. However, the angular 
momentum quantum numbers J and K typically have well-defined mean values that are 
useful for interpretation. 

Figure El shows surface plots of the mean values (J) and {K) for the FL potential surface 
(note that the axes are rotated relative to Fig. 12}. Near the minima of the potential wells. 
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M 


V 


Energy (K) 








0.0282 





1 


0.00550 





2 


0.000455 


1 





0.00545 



TABLE I: Binding energies in Kelvin of FL states. Each state is identified by its rotation M about 
the electric field axis, and by a vibrational quantum number v. These energies refer to states even 
under the reflection 9 ^ it — 9. Additional states, odd under this symmetry, are separated in 
energy by less than several fiK from the ones listed. 




FIG. 3: Average values of J and K for the Ai = potential in Fig. 2. Note that the orientation 
of the axes is different from that in Figs 2 and 4. 

we find that (J) ~ (K) ^ 2. These values characterize the FL states at large separation 
R. Ref. 0] presents a simple and useful semiclassical picture of the dipole's orientation 
in the OH molecule. In this model the dipole precesses around the molecule's total angular 
momentum j, and on average points along j when uj > 0, and against j when uj < 0. Thus 
when ( J) ^ (K) ^ 2 and ooi 002, as is the case here, the dipole moments are both aligned on 
average in the same direction, roughly along the intermolecular axis, and hence attract one 
another. 

At smaller values of R, (J) remains nearly equal to 2, but (K) drops all the way to 
0. This reflects the influence of the avoided crossings in the surfaces. Again invoking a 
semicalssical picture, (J) = 2, (K) = implies that the dipole moments are now aligned 
roughly perpendicular to the intermolecular axis, in a side-by-side orientation where they 
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repel one another. This is the reason the FL state is stable against collapse to smaller R. 

The avoided crossings that allow FL states to be supported have their origin in the fact the 
the Stark interaction is diagonal in the laboratory frame (defined by the field axis), whereas 
the dipolar interaction is diagonal in the rotating frame (defined by the intermolecular axis). 
Competition between these two symmetries generate the avoided crossings. However, in the 
limit where 6^0 the two axes coincide and both interactions become diagonal in K. In 
this case the avoided crossings become diabatic crossings, and there is a conical intersection 
in the surfaces. Our description in terms of adiabatic surfaces is, therefore, incomplete. It 
is however useful, as we will see in the next section. There may be interesting information 
on geometrical phases inherent in the FL states; this will be a topic of future study. 



B. Bound states 



To complete a description of the FL states we must understand their motion in R and 6. 
Each bound state is nearly doubly-degenerate with respect to reflection in the 6 = n /2 plane. 
In Figure m we present wave function plots of those bound states that have even reflection 
symmetry, corresponding to the bound states listed in Table I. In this figure, (a-c) refer to 
the Ai = case, and (d) to the A4 = 1 case, for = 0, it is immediately evident that 
these states exhibit zero-point motion in the 6 direction, and that excitations are primarily 
in the R direction. We therefore label the FL states with a vibrational quantum number v. 
For Ai = 1, a nodal line appears along the ^ = direction, owing to the centrifugal energy 
that forces the molecules away from the electric field axis. 

In realistic laboratory circumstances, the FL states are quasi-stable, being subject to 



dissociation into free molecules in lower-energy internal states [17[. Nevertheless, the adia- 
batic bound states we have identified here correspond to real features of these dissociating 
states. To show this, we have carried out a complete coupled- channel scattering calculation 



in a laboratory- frame representation, similar to that in Ref. la], but without including 
hyperfine structure. We have included partial waves up to L = 16 to ensure convergence at 
the several percent level in scattering observables. 

We compute the time delay for the scattering process, defined as 

r^2n^, (21) 
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FIG. 4: Wave functions of FL states, for the potential surfaces shown in Fig. 2. For = 0, there 
are three vibrational states, u = (a), u = 1 (b), and v = 2 (c). For = 1, there is a single state, 
with V = (d). 




Binding Energy (K) 

FIG. 5: Time delay as defined in Eqn. (22) for the = (a) and M = 1 (b) FL states. 
Resonance peaks appear at cliaracteristic energies tliat correlate well with the binding energies as 
determined from the FL adiabatic surfaces (vertical lines). 

where 6 is the eigenphase sum, i.e., the sum of the inverse tangents of the eigenvalues of the 
scattering K-matrix. This quantity, plotted in Figure exhibits peaks at energies where 
resonances occur. There is also a significant background component, arising from threshold 
effects, but the peaks are nevertheless visible. Also indicated are the binding energies of the 
FL states as given in Table I. The good agreement between the two calculations suggests 
that when FL resonances are observed in experiments, the nature of the resonant states will 
be well-approximated by the wave functions determined above. 

In general, resonances in ultracold polar molecule scattering will come in three varieties. 
The "true" field-linked states, like the ones that we describe here, are largely independent of 
physics at small values of R. We can verify this assertion by changing the small- boundary 
conditions in our multichannel scattering calculation. The positions of the FL resonances 
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do not depend at all on these boundary conditions. However, their lifetimes can fluctuate 
within a factor of ~ 2, since the continuum states into which they can decay do depend on 
short-range physics. 

A second type of FL state appears to have components at both large and small R. 
Examples of these are found for states lying below the middle threshold in Fig. ^ We flnd 
that their positions are relatively insensitive to the short-range boundary conditions, but 
that their lifetimes vary wildly. We refer to these as "quasi-FL" states. Finally a third 
category of resonance is strongly sensitive to initial conditions, both in position and width. 
These are resonant states of the short-range interaction, which are expected to be numerous 
in low-energy molecular collisions 0, Is^l . 

IV. OUTLOOK 

We have left out many details of molecular structure and interactions, in order to em- 
phasize the basic structure of the fleld-linked states. This structure is remarkably simple, 
and consists primarily of a pair of molecules in relative vibrational motion along an axis 
that nearly coincides with the direction of the electric fleld. The number of FL states is not 
large, since the forces holding them together are necessarily weak. 

Signiflcantly, to adapt this simple picture to a particular molecular species requires only 
a detailed knowledge of the structure of each molecule separately, plus some information on 
long-range parameters such as dispersion coefficients. In other words, realistic modeling of 
experimentally probed FL states can probably be achieved using currently existing informa- 
tion. This is in stark contrast with molecular collisions involving close contact between the 
molecules, in which case existing potential energy surfaces are likely to be inadequate for to 
describe collisions at ultralow temperatures. 
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APPENDIX A: SYMMETRIZED WAVE FUNCTIONS 



To incorporate the effects of symmetrization under the exchange (P12) and parity (/) 



operations, we follow the treatment of Alexander and DePristo [27[. To this end it is 
convenient to relate the Euler angles of each molecule to the electric field axis rather than 
the intermolecular axis; these Euler angles are denoted e^. The symmetry operations then 
perform the following functions: 



Pi 



12 



R 



R 



I: R^-R 

e[ /(ef ) 

The last two lines imply that / acts on each molecule by inverting the molecule's coordinates 
through its own center of mass. 

The effect of particle exchange on the internal coordinates is determined by making the 
explicit rotation to the lab frame: 

A2(ei,e2|(l,2)Ji^) 

= A2 E (ef , e^l (1, 2) Jmi2)I^4„^(</>, 0) 

mi2 

= ^(_l)W(ef,e^|(2,l)Jmi2)(-l)^/^4„_^(0,^,O) 

mi2 

= (-lf^-(ei,e2|(2,l)J-ir). 



Here we have used the reflection symmetry of the Wigner D functions, 

Dij,{T: + 0, vr - ^, 0) = {-iyDi_K{ct^, 0) 

and the usual exchange symmetry of the Clebsch-Gordan coefficients. Similarly the relative 
wave functions transform as 
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exp(zX0)F(^2)^^(P,^) 
(-l)^exp(zA^0)F(-t;'2)j;.(^,^ 
17 



An appropriately symmetrized basis for exchange is therefore given by Eqn. where 



^(1 



(1,2)JK ~ 2 

X 



X 



^^2(1 + 512^ 
l,2)JK)+x(-lf^|(2,l)J-ir) 



with sx = ±1 for bosons/fermions. 

These basis functions can in turn be assembled into parity eigenf unctions. Note that / 

has the same effect on the relative coordinates as does Pu, so that g'^^'^F^2)jk already 

a parity eigenstate with eigenvalue s. Denoting the parity of the total wave function by e, 

the parity of the relative wavefunctions should be p = es, or p = q{—l)^ in terms of our 

quantum number q defined in Eqn.®. This definition seems (and is) completely arbitrary; 

it is justified by explicitly working out the matrix elements for the Stark and dipole-dipole 

interactions, and finding that both conserve the value of q. 

The influence of I on each molecule is to reverse its direction of rotation about its own 

I — I 

axis, and to introduce a phase |2Si] : 



/(e^|j,m,^) = (-l)^-"(e^|jm-^). 

Because the phase factor is the same for each molecule, the action of / on the molecule pair 
is, by arguments similar to those above, 

I{e^e2\{l,2)JK) = {-ly {e.e^H-l, -2)J - K), 

where the notation (1,2) (—1,-2) implies (ji, cji, j2, ^^2) (ji, — ^^i, j2, — ^^2)- The sym- 
metrized internal basis function is then 



:i,2)ji^).,, 
1 

72 



1(1, 2) JK), + g(-l)^+^'|(-l, -2) J - K), 



APPENDIX B: CONSERVATION OF q 



It is straightforward (if somewhat tedious) to write the symmetrized matrix elements for 
different contributions to the Hamiltonian, in terms of the unsymmetrized basis. We present 
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here some of the key results, which rely mostly on the symmetry properties of the angular 
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momentum recoupling coefficients, as described in Brink and Satchler 

Dipolar interaction. In Eqn. (I15|) . the 9-j symbol must be invariant under exchanging its 
second and third rows, yet this operation introduces a phase shift ( — 1)"'+^+"''. Therefore, 
we must have J+ J' = even, and the matrix element (|15|) is invariant under the substitution 
(1, 2) (2,1), K ^ —K. In the symmetrized basis it reads 

^{{l,2)JK\V^^\{l\2')J'K'),= 
11 + xx'\ 1 



X 



2 ) \J{1 + 5i25kq){1+5v2'5k'q) 
{{l,2)JK\V,,\{l\2')J'K') 

+x'{-lf^{{l,2)JK\V,,\{2\l')f-K') 

Thus the exchange quantum number x is explicitly conserved. Similarly, the matrix elements 
are invariant under simultaneously reversing the signs of all cu's and /T, whereby 

,,,((l,2)Jir|V^^|(l',2')J'i^%V = 

2 J 

{{l,2)JK\V^^\{l\ 2')J'K')^, 

-X 

+g'(-l)''' 1)JK\V,,\{-V, -2')J' - K'),, 

Because K = K' for the dipolar interaction, this implies in turn that q is conserved. The 
matrix derivation of symmetrized matrix elements for the lamb da- doubling is exactly the 
same, and this interaction also conserves q. 

Stark interaction. Symmetrized matrix elements of the Stark Hamiltonian (fTT|) are 
slightly more complicated, since reversing the sign of K also affects the Wigner (i-function. 
Exploiting symmetries of the d functions yields 



,{{l,2)JK\H\l\2')J'K') 



XX 



\K+K' 



1 



{{l,2)JK\H'\{l',2')J'K')d'M 



+x\-lYHil,2)JK\H^\i2',l')J' - K')dl^,^M 
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In general, neither x, nor K, nor the product x{—l)^, is conserved by this part of the 
Hamiltonian. However, the matrix elements in the basis (|7j) become 

,,((l,2)Jir|ff^|(l',2')J'i^%W 
{il,2)JK\H'\{l',2')J'K'),> 



X 



+q'i~lf .'((1, 2)JK\H'\i-l', -2')J' - K'), 

illustrating the conservation of q. 

Centrifugal energy. Symmetrized over x, the centrifugal energy reads 



< {l,2)JK\V^'"^'\{l',2')J'K' 



SjT' X 



1 + xx'^ 



2mi?2"^'- V 2 ; V^4(l + 5i,25ir,o)(l + 5r,2'5A",o)/ 



+x'(Ao{K)5.K,K' + A^i{K)6-K,K'+i + A+i{K)6 -K,K'-1 )5l,2'52,l 



From this point, translation into the x, q symmetrized basis is trivial. In general, q is not 
conserved by the Coriolis terms that change K, but in the present treatment these terms 
are ignored. 
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